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An  accurate  assessment  of  liquid-phase  mass  transport  resistances  is  necessary  for  understanding  and 
optimizing  battery  performance  using  mathematical  models.  This  work  combines  modeling  and  experi¬ 
ments  to  quantify  tortuosity  in  electrolyte-filled  porous  battery  structures  (separator  and  active-material 
him).  Tortuosities  of  separators  were  measured  by  two  methods,  AC  impedance  and  polarization- 
interrupt,  which  produced  consistent  results.  We  measured  an  apparent  interfacial  resistance  at  the 
lithium  metal  electrodes  that  contributed  to  both  ohmic  and  diffusional  resistance  of  the  cell.  The 
polarization-interrupt  experiment  was  used  similarly  to  measure  effective  electrolyte  transport  in  porous 
films  of  cathode  materials,  particularly  films  containing  LiFeP04.  An  empirical  relationship  between 
porosity  and  the  tortuosity  of  the  porous  structures  was  developed.  Our  results  demonstrate  that  the 
tortuosity-dependent  mass  transport  resistance  in  porous  separators  and  electrodes  is  significantly  higher 
than  that  predicted  by  the  oft-used  Bruggeman  relationship. 

©  2008  Elsevier  B.V.  Ah  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  widely  used  in  consumer  electronic 
devices.  New  applications  like  hybrid  electric  vehicles  (HEVs)  and 
plug-in-hybrid  electric  vehicles  (PHEVs)  are  potential  candidates 
for  using  Li-ion  batteries.  To  achieve  required  energy  and  power 
densities,  it  is  necessary  to  optimize  batteries  by  understand¬ 
ing,  quantifying,  and  mitigating  the  dominant  resistances.  One 
approach  to  increased  understanding  has  been  the  use  of  macro- 
homogeneous  mathematical  models  [1-9].  The  effectiveness  of 
such  models  depends  to  a  large  extent  on  the  accuracy  of  the  phys¬ 
ical  parameters  used.  The  purpose  of  this  work  is  to  examine  and 
increase  the  accuracy  of  the  parameters  that  describe  liquid-phase 
mass  transport  in  electrodes  and  separators  used  in  Li-ion  batteries. 

Porous  electrodes  and  separators  are  multi-phase  structures 
that  include  a  network  of  interconnected  and  irregular  pores  or 
channels.  It  is  difficult  and  time  consuming  to  model  these  compli¬ 
cated  pore  networks  in  full  three-dimensional  realism.  Instead,  for 
macro-homogeneous  models,  the  complexity  of  mass  transport  in 
the  pore  network  is  reduced  to  only  one  effective  geometric  param¬ 
eter,  tortuosity  r.  The  tortuosity  of  a  porous  network  is  used  to 
obtain  effective  transport  properties.  It  is  assumed  that  the  tor¬ 
tuosity  affects  transport  properties,  liquid-phase  diffusivity  and 
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conductivity,  with  a  similar  functionality  [10].  In  general, 
kc 

Keff  =  ~  (1) 

or 


Here  /ceff  and  /c  are  the  effective  and  intrinsic  conductivities,  respec¬ 
tively,  and  Deff  and  D  are  the  effective  and  intrinsic  diffusivities, 
respectively,  of  the  conductive  phase  (in  the  present  case  the  liquid 
electrolyte).  Eqs.  (1 )  and  (2)  serve  as  the  definition  of  r  used  in  the 
present  work.  When  comparing  r  values  from  different  sources,  one 
must  be  cognizant  of  the  definitions  used.  For  example,  our  tortu¬ 
osity  is  the  square  of  the  one  used  by  Abraham  [11  ].  Alternatively, 
some  workers  choose  to  describe  effective  conductivity  in  terms  of 
a  MacMullin  number  [12],  which  in  this  context  is  Nm  =  r/c. 

The  effective  transport  properties  are  corrected  by  e,  the  volume 
fraction  of  the  conductive  phase  (in  the  present  case  €  is  the  porosity 
or  void  fraction  of  the  solid).  This  means  that  any  salt  molar  fluxes 
(from  Deff)  or  electrical  current  densities  (from  /ceff)  in  the  model  are 
based  on  superficial  electrode  and  separator  areas.  If  the  conductive 
pathways  are  composed  of  straight  channels  of  uniform  cross  sec¬ 
tion  that  are  parallel  to  the  transport  direction,  then  r  =  1  in  Eqs. 
(1)  and  (2).  Otherwise,  for  tortuous  pathways  r  >  1  and  effective 
transport  properties  are  always  lower  than  the  intrinsic  transport 
properties  of  the  conductive  phase. 

Of  course,  r  depends  in  a  complicated  way  on  the  structure  of  the 
porous  network,  including  the  nature  and  particle-size  distribution 
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of  each  of  the  materials  used  in  the  porous  electrode.  Due  to  these 
complexities,  researchers  have  used  r  as  an  adjustable  parameter 
or  perhaps  as  a  function  of  e  [10,13].  The  Bruggeman  relationship 
[11,14-17]  has  often  been  used  to  correlate  and  predict  effective 
conductivities: 

Keff  =  /ee“,  (3) 

or  equivalently, 

r  =  e1-“,  (4) 

where  a  is  the  Bruggeman  exponent.  For  a  system  composed  of  a 
continuous  conductive  phase  (e.g.  liquid  electrolyte)  mixed  with 
insulating  spherical  particles  of  uniform  size  (i.e.  monodisperse)  it 
has  been  shown  empirically  [16,17]  that  o'  ^  1.5.  This  value  of  a  is 
often  used  in  battery  simulations. 

There  is  almost  a  complete  lack  of  direct  experimental  data  for 
tortuosity  or  liquid-phase  transport  resistance  in  porous  electrodes. 
Similarly,  few  have  questioned  the  accuracy  of  the  Bruggeman 
equation  (Eq.  (3))  with  a  =  1.5.  A  few  researchers  [10-12]  have 
used  larger  values  of  the  Bruggeman  exponent.  Patel  et  al. 
[12]  performed  numerical  simulations  of  porous  networks  with 
different  morphologies  and  calculated  the  tortuosity  for  these 
networks  in  terms  of  a  Bruggeman  exponent.  They  also  per¬ 
formed  liquid-phase  conductivity  measurements  for  different  types 
of  battery  separators  and  calculated  corresponding  Bruggeman 
exponents. 

Prior  work  by  our  group  [18]  found  that  liquid-phase  transport 
resistance  in  LiCo02  cathodes  has  a  great  effect  on  cell  performance 
and  that  the  usual  Bruggeman  relation  with  a  =  1.5  significantly 
underpredicts  mass-transport  resistances.  While  one  can  regress 
needed  physical  parameters  by  comparing  cell  performance  to  a 
full-cell-sandwich  model  (as  we  did),  the  presence  of  multiple 
unknown  resistances  (including  kinetic  and  electronic  transport) 
in  the  model  makes  accurate  determination  of  one  particular  resis¬ 
tance  difficult.  Therefore,  it  is  of  great  interest  to  develop  an 
experimental  approach  to  quantify  directly  the  tortuosity  or  the 
liquid-phase  transport  resistance  for  porous  electrode  films  and 
separators.  We  are  not  aware  of  prior  experimental  evaluation  of 
the  tortuosity  of  porous  electrodes  independent  of  other  processes 
in  the  cell. 

In  this  work  we  present  techniques  to  determine  the  tortuosities 
of  different  Li-ion  active-material  films  (LiFeP04  and  LiCo02)  and 
one  type  of  battery  separator  (Celgard  2400).  Experiments  included 
AC  impedance  to  obtain  /ce ff  and  a  polarization-interrupt  method  to 
obtain  Deff,  each  performed  on  cells  with  a  symmetric  geometry 
and  using  an  electrolyte  with  known  intrinsic  transport  proper¬ 
ties  k  and  D.  The  computer  model  described  below  was  developed 
for  interpreting  the  experiments  and  is  essential  for  the  success 
of  the  method.  The  use  of  two  complementary  experiments,  one 
for  Ke ff  and  one  for  Deff,  and  of  different  thicknesses  of  materials 
serve  to  validate  the  results.  Finally,  we  present  a  summary  plot  of 
the  results  showing  the  effect  of  porosity  on  the  tortuosity  of  the 
active-material  films. 


2.  Correlating  tortuosity 

As  indicated  by  the  above  correlations  (Eqs.  (1)  and  (2))  the 
porosity  e  plays  an  important  role  in  the  electrochemical  perfor¬ 
mance  of  an  electrode.  Lower  porosity  corresponds  to  a  greater 
fraction  of  solid  particles  including  carbon,  and  hence  a  greater 
number  of  possible  electronic  pathways.  This  means  that  lower 
porosity  is  desired  for  lower  electronic  resistance  [19].  On  the  other 
hand,  higher  porosity  or  liquid  volume  fraction  is  desired  for  lower 
liquid-phase  transport  resistance.  Therefore,  to  optimize  battery 
for  a  particular  application,  it  is  important  to  know  the  effect  of 
porosity  on  the  liquid-phase  resistance  in  the  electrode. 

In  this  work,  we  use  a  generalized  form  of  the  Bruggeman  rela¬ 
tion  [14]  to  correlate  tortuosity: 

T  =  y€1~a,  (5) 

where  y  and  a  are  the  constants  that  could  depend  on  the  morphol¬ 
ogy,  porosity,  material,  and  particle-size  distribution  of  a  porous 
composite.  Table  1  summarizes  some  prior  published  experimen¬ 
tal  work  in  terms  of  Eq.  (5).  As  shown,  most  of  the  prior  work 
involves  either  setting  y  =  1  and  adjusting  a,  or  setting  a  =  1 .5  and 
adjusting  y,  in  order  to  obtain  the  needed  r  value.  Without  multiple 
experiments  at  different  e  values,  one  cannot  determine  the  proper 
relationship  between  r  and  e,  if  any. 

Patel  et  al.  [12],  Djian  et  al.  [20],  and  Abraham  [11]  calculated 
tortuosities  for  different  types  of  separators.  For  Celgard  2400  sep¬ 
arators,  these  groups  reported  different  values  for  the  tortuosities. 
However,  the  Celgard  2400  samples  used  by  these  groups  appear 
to  have  different  porosities.  The  separator  used  in  this  work  has  the 
same  porosity  as  that  used  by  Patel  et  al. 

Table  1  also  compares  tortuosities  calculated  for  porous  elec¬ 
trodes  and  for  a  packed  bed.  The  Bruggeman  approximation  was 
reconfirmed  by  Patel  et  al.  [12]  and  Pavlin  et  al.  [21]  for  conduc¬ 
tivity  in  a  continuous  phase  mixed  with  uniform-size  spherical 
particles.  However,  from  numerical  simulations  for  different  par¬ 
ticle  shapes,  Patel  et  al.  [12]  found  that  the  tortuosity  was  higher 
than  that  predicted  by  the  Bruggeman  relation.  Also,  Doyle  et  al. 
[10]  for  LiMn204  electrodes  and  Stephenson  et  al.  [18]  for  LiCo02 
electrodes  found  that  the  tortuosity  for  the  porous  electrodes  is 
significantly  higher  than  the  traditional  Bruggeman  relationship. 
However,  these  two  works  used  discharge  curves  to  fit  r  as  an 
adjustable  parameter.  The  cell  examined  by  Doyle  et  al.  contained 
a  gelled  electrolyte,  consequently  the  reported  tortuosity  accounts 
not  only  for  the  pores  between  active-material  particles,  but  also 
for  the  transport  retardation  effect  from  intimate  mixing  of  poly¬ 
mer  and  liquid  electrolyte.  As  researchers  have  used  tortuosity  of 
cathode-material  films  as  an  adjustable  parameter  in  their  full  bat¬ 
tery  simulations,  the  reported  tortuosity  values  may  inadvertently 
incorporate  the  effects  of  other  processes  (kinetics,  transport,  and 
thermodynamics)  not  accurately  captured  by  the  model.  To  accu¬ 
rately  quantify  the  liquid-phase  transport  resistance  or  tortuosity, 
it  is  imperative  to  isolate  the  liquid-phase  transport  processes  from 
the  other  processes  in  the  electrode.  This  work  presents  experi¬ 
mental  and  modeling  techniques  that  are  developed  to  separate  as 
much  as  possible  the  liquid-phase  transport  resistances  from  other 
processes. 


Table  1 

Reported  tortuosities  in  the  form  z  =  ye1_“  for  selected  porous  structures. 


Author 

y 

a 

r 

6 

Porous  structure 

MacMullin  and  Muccini  [17] 

i 

1.6-1. 7 

4-5 

0.1 

Bed  of  spherical  sandstone  particles 

Abraham  [11] 

7.3 

1 

7.3 

0.32 

Separator  (Celgard  2400) 

Djian  et  al.  [20] 

1 

2.46 

5.3 

0.32 

Separator  (Celgard  2400) 

Patel  et  al.  [12] 

1 

2.8 

6.0 

0.37 

Separator  (Celgard  2400) 

Doyle  et  al.  [10] 

1 

3.3 

2.9 

0.63 

LiMn204  electrode 

Stephenson  et  al.  [18] 

18 

1.5 

21-27 

0.43-0.71 

UC0O2  electrode 
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3.  Experimental  procedure 

This  section  describes  the  experimental  techniques  used  to  fab¬ 
ricate  cells  and  quantify  liquid-phase  transport  processes  for  a 
separator  (Celgard  2400)  and  for  free-standing  porous  cathode 
films. 

3  A.  Cell  fabrication 

To  quantify  the  tortuosity  of  the  porous  structures,  a  symmet¬ 
ric  cell  geometry  was  used  that  is  consistent  with  the  method  of 
restricted  diffusion.  Restricted-diffusion  experiments  involve  one¬ 
dimensional  diffusion  in  cells  closed  at  the  ends  and  in  which 
concentration  variations  take  place  over  the  entire  length  between 
the  ends  [22,23].  Fig.  1  is  a  schematic  of  the  cell  geometry  used  to 
determine  the  tortuosity  of  the  separator.  The  cells  were  assem¬ 
bled  in  an  Ar-atmosphere  glove  box  (moisture  content  0.9  ppm 
and  oxygen  content  <  0.25  ppm,  VAC,  Hawthorne,  CA).  Lithium 
metal  foil  (Alfa  Aesar)  of  area  4  cm2  on  a  copper  current  collec¬ 
tor  was  used  for  each  electrode.  The  two  lithium  electrodes  were 
separated  by  25-p.m-thick  polypropylene  separator  layers  (Cel¬ 
gard  2400,  Celgard  LLC,  Charlotte,  NC).  Metalized  polymer  film 
(Class  PPD  Shield  Pack,  Inc.,  Los  Angeles,  CA)  was  used  to  enclose 
the  cells.  Electrolyte  composed  of  1  M  LiPF6  in  a  1:1  (w:w)  mix¬ 
ture  of  ethylene  carbonate  and  diethyl  carbonate  (LithDyne  Elyte, 
LithChem  International,  Anaheim,  CA)  was  added  to  the  pouch 
before  sealing  with  an  electric  heat  sealer  (Impulse  Sealer  Tish 
200,  Electronic  Heating  Equipment  Co.).  The  cells  were  tested 
under  external  compression  of  68kPa  above  atmospheric  pres¬ 
sure. 

A  geometry  similar  to  Fig.  1  was  used  to  quantify  the  trans¬ 
port  properties  for  active-material  films  as  shown  in  Fig.  2— the 
main  difference  is  the  presence  of  a  free-standing  active-material 
film  placed  between  the  two  lithium  electrodes.  Separator  layers 
electronically  isolated  the  active-material  film  from  the  lithium 
electrodes  as  shown  in  Fig.  2.  The  active-material  films  were  fabri¬ 
cated  with  one  of  two  different  active  materials  (LiFeP04  or  LiCo02) 
and  delaminated  from  their  current  collectors  as  described  below. 
In  this  paper  by  convention  we  refer  to  the  free-standing  active- 
material  films  or  cathode  films  based  on  their  composition  and 
structure,  not  their  function;  the  true  electrodes  in  our  experimen¬ 
tal  cells  were  always  a  pair  of  lithium  foils. 

Each  active-material  film  contained  (on  a  dry  basis)  84  wt% 
active  material,  either  LiFeP04  (Hydro  Qubec,  Montral,  Canada) 


Fig.  1.  Schematic  of  a  symmetric  cell  geometry  for  testing  separators.  Lithium  foils 
are  used  as  anode  and  cathode,  between  which  are  placed  one  or  more  separator 
layers  filled  with  the  electrolyte. 
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Fig.  2.  Schematic  of  a  symmetric  cell  geometry  for  testing  porous-electrode  films. 
The  active-material  film  is  delaminated  from  its  current  collector  and  placed 
between  separators,  which  in  turn  are  placed  between  lithium  foil  electrodes. 

or  LiCo02  (Alfa  Aesar),  along  with  4  wt%  carbon  black  (diameter: 
20-30  nm,  XC-72R,  Cabot  Carbon  Ltd.),  4wt%  graphite  (diame¬ 
ter:  5. 5-7.5  |jim,  SFG-6,  TIMCAL  Group),  and  8  wt%  polyvinylidene 
fluoride  copolymer  (PVdF,  741  Kynar  Corp.).  In  the  case  of  carbon- 
fiber-containing  LiFeP04-material  films,  the  graphite  and  carbon 
black  components  (8wt%  total)  were  replaced  by  carbon  fibers 
(diameter:  100-200  nm,  length:  30-100  [xm,  PR-19,  Pyrograf  Prod¬ 
ucts  Inc.).  To  fabricate  the  active-material  film,  PVdF  was  dissolved 
in  NMP  solvent  (l-methyl-2pyrrolidone  fluoride)  with  mild  heat¬ 
ing.  The  active-material  and  carbon  components  were  dry-mixed 
with  a  spatula  followed  by  mixing  with  an  ultrasonic  homogenizer. 
This  solid  mixture  was  then  mixed  with  the  solution  of  PVdF  and 
NMP  with  a  spatula  followed  by  mixing  with  the  ultrasonic  homog¬ 
enizer  to  form  a  composite  slurry.  Then,  this  composite  slurry  was 
spread  onto  an  aluminum  current  collector  with  a  doctor  blade 
(Byk-Gardner,  USA)  to  a  desired  wet  thickness.  To  facilitate  later 
delamination  the  slurry  was  applied  to  the  shiny  (thicker  oxide) 
side  of  the  aluminum  foil  without  the  benefit  of  an  adhesion  layer. 
The  composite  cast  was  dried  overnight  at  120°  C  under  vacuum 
(15  mmHg). 

Individual  active-material  films  were  cut  from  the  dried  com¬ 
posite  cast  using  a  template  of  10  cm2  area.  The  area  of  the  active 
material  films  was  larger  than  that  of  the  lithium  electrodes  to 
make  sure  that  transport  of  all  lithium  ions  between  electrodes 
took  place  only  through  the  active  material  film.  The  active-material 
films  were  compressed  to  different  thicknesses  using  a  calendering 
machine  (Pepetools,  OK).  From  weight  and  thickness  measure¬ 
ments,  we  calculated  the  porosities  of  the  active-material  films. 
After  the  calendering,  the  aluminum  current  collector  was  care¬ 
fully  removed  from  these  active-material  films  using  a  0.009-inch 
single-edge  razor  blade  as  shown  in  Fig.  3,  with  care  taken  to  min¬ 
imize  bending  of  the  films.  The  delaminated  active-material  film 
was  then  assembled  into  a  pouch  cell  as  previously  discussed.  Once 
the  cell  was  fabricated,  it  was  tested  using  the  procedures  described 
below. 

3.2.  Cell  testing 

The  tortuosity  of  the  separator  or  active-material  films  was 
obtained  from  two  methods:  polarization-interrupt  (or  restricted 
diffusion)  and  AC-impedance.  The  polarization-interrupt  test  was 
performed  on  assembled  cells  to  obtain  the  effective  diffusivity 
of  the  electrolyte.  Similarly,  AC  impedance  was  performed  on  the 
assembled  cells  to  calculate  the  effective  conductivity  of  the  elec¬ 
trolyte.  In  each  case  the  effective  diffusivity  or  conductivity  was 
used  to  calculate  the  tortuosity  of  the  porous  structure  using  Eqs. 
(1)  or  (2). 
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Fig.  3.  Image  showing  procedure  to  remove  aluminum  current  collector  (curled  top 
layer)  from  active-material  film  (bottom  black  layer)  using  a  razor  blade. 

3.2.1  Polarization-interrupt  experiment 

The  assembled  cells  (separator  and  active-material  film  cells) 
were  tested  using  a  battery  tester  (Maccor  4300,  Maccor  Inc.).  Sev¬ 
eral  conditioning  cycles  were  performed  on  a  cell  to  form  a  stable 
SEI  (solid  electrolyte  interface)  on  the  lithium  electrodes.  Each  cycle 
consisted  of  passing  a  constant  current  (0.5  mA  cm-2)  for  10  min, 
followed  by  a  3-min  rest  period,  then  passing  current  in  the  reverse 
direction  under  the  same  conditions.  After  the  conditioning  cycles, 
a  constant  current  was  applied  for  2  min  to  set  up  a  concentration 
gradient  in  the  cell.  Then,  the  concentration  gradient  was  allowed  to 
relax  by  interrupting  the  current  until  the  cell  potential  approached 
zero.  Then  the  current  direction  was  reversed  for  2  min  followed 
by  interruption  of  the  current  until  cell  potential  approached  zero. 
These  steps  were  repeated  for  different  current  densities.  The  relax¬ 
ation  of  the  concentration  gradient  was  observed  in  terms  of  a 
cell-potential  decay  with  the  time  (Figs.  6  and  8).  As  discussed  in 
more  detail  below,  the  model-experiment  comparison  (discussed 
in  Section  5 )  of  this  cell-potential  decay  behavior  was  used  to  obtain 
quantitative  information  about  the  tortuosities  of  the  porous  struc¬ 
tures. 

3.2.2.  AC  impedance  experiment 

After  the  polarization-interrupt  experiment,  a  current- 
controlled  AC  impedance  test  was  performed  using  a  frequency 
response  analyzer  (Gamry  Instruments)  on  the  separator  cells  to 
measure  the  total  resistance  between  the  two  electrodes,  R^.  The 
frequency  of  the  AC  signal  was  varied  from  100  Hz  to  100  kHz  at 
an  amplitude  of  0.025  mA  cm-2.  was  taken  from  the  intercept 
with  the  real  axis  at  high  frequency.  Fig.  4  shows  a  representative 
Nyquist  plot  obtained  from  the  AC  impedance  experiment  for  one 
of  the  cells.  R^  is  then  used  in  conjunction  with  a  model  (see 
below)  to  obtain  tortuosity  for  the  separator. 

4.  Model  development 

In  order  to  interpret  each  type  of  experiment  and  obtain  a  tor¬ 
tuosity  value,  the  following  models  were  used. 

4.1  Polarization-interrupt  experiment 

The  equations  describing  transport,  kinetic,  and  thermody¬ 
namic  processes  in  porous  electrodes  have  been  outlined  previously 
[24-26].  In  this  work,  during  the  polarization  of  a  cell,  lithium  ions 
are  generated  at  the  anode  surface,  diffuse  and  migrate  through 
the  porous  regions  (separators  and  active-material  film),  and  are 
deposited  on  the  cathode.  In  order  to  model  this  system,  the  gov¬ 


Fig.  4.  Nyquist  plot  obtained  from  AC  impedance  experiment  for  separator  cell 
shown  in  Fig.  1.  The  intercept  with  the  real  axis  at  high  frequency  (R^)  is  used 
to  find  the  effective  ionic  conductivity  in  the  separator  region. 

erning  equations  for  the  processes  in  the  porous  regions  were 
implemented  and  solved  using  a  finite  element  package  (COMSOL 
Multiphysics  3.4)  as  described  previously  by  Stephenson  et  al.  [18]. 
However,  because  the  cell  geometry  was  not  the  same  as  that  mod¬ 
eled  by  Stephenson  et  al.,  modifications  to  model  geometry  and  the 
governing  equations  were  made. 

The  separator  and  active-material  film  regions  were  modeled 
with  a  ID  restricted-diffusion  geometry.  The  lithium  metal  anode 
and  cathode  (see  Fig.  2)  appear  as  boundary  conditions  in  the  model 
(x  =  0,  L ).  A  brief  summary  of  the  governing  equations  and  bound¬ 
ary  conditions  used  in  the  model  follows. 

4  A  A.  Charge  transfer  through  the  liquid  phase 

For  the  electrolyte  phase,  the  governing  equation  for  the  porous 
regions  is 


where  i  is  the  superficial  current  density  in  the  liquid  phase  and  x 
is  the  position  in  the  direction  normal  to  the  electrode  surfaces. 

We  assume  a  binary  electrolyte,  no  convection,  and  concen¬ 
trated  solutions  with  a  constant  cation  transference  number  t°. 
Given  this,  the  current  in  the  electrolyte  can  be  expressed  in  a 
particular  porous  region  by 

3 4>  2KeffRT  ,  o x  f  d\nf±\  3 Inc 

1  =  -*eff3^  +  — f^W+H1+  dHFj^T’  (7) 

where  0  is  the  electric  potential,  R  is  the  ideal  gas  constant,  T  is 
the  temperature,  F  is  the  Faraday’s  constant,  and  c  is  the  elec¬ 
trolyte  concentration.  The  appropriate  effective  conductivity,  xef f, 
for  each  porous  layer  was  used.  The  activity  coefficient  correction, 
1  +  d  In  f±/d  In  c,  in  this  work  was  taken  to  be  one. 

In  a  recent  paper,  Stewart  and  Newman  [27]  measured  the 
activity  coefficient  correction  for  LiPF6  in  three  different  carbon¬ 
ate  solvents,  though  not  the  mixed  solvent  used  in  this  work.  In 
each  case  the  values  increased  with  concentration,  reaching  a  value 
around  2.5  in  the  case  of  1  mol  dm-3  LiPF6  in  1:1  (w:w)  EC/EMC. 
The  polarization  experiments  in  this  work  involved  relatively 
small  concentration  variations  about  the  average  concentration 
of  1  mol  dm-3.  Inclusion  of  a  concentration-dependent  correction 
factor  in  Eq.  (7)  therefore  does  not  substantially  influence  the  mea¬ 
sured  diffusive  time  constant  of  the  relaxation  period  following 
polarization,  and  hence  the  calculated  tortuosity  results.  All  the 
potentials  were  taken  relative  to  the  lithium  anode.  The  boundary 
condition  at  the  anode  surface  (x  =  0)  is  given  by  a  Butler-Volmer 
expression  (Eq.  (6)  in  Ref.  [18]).  The  exchange  current  density  for 
the  lithium  metal  electrodes  was  obtained  by  the  AC-impedance 
method  and  was  set  to  io  =  0.41  mAcm_2at  c  =  1  mol  dm-3.  The 
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boundary  condition  at  the  cathode  side  (x  =  I)  is  i  =  /,  where  I 
is  the  fixed  current  density  being  passed  through  the  cell.  The 
liquid-phase  potential  and  superficial  current  are  continuous  at  the 
porous-layer  interfaces.  A  Butler-Volmer  expression  at  the  cathode 
electrode  surface  provides  for  the  potential  of  the  cathode  relative 
to  the  anode.  In  effect  the  model  sums  the  two  kinetic  overpoten¬ 
tials  and  the  ohmic  drop  and  concentration  overpotential  of  each 
layer  to  obtain  the  cell  potential. 

4.1.2.  Diffusion  in  the  liquid  phase 

The  continuity  equation  for  mass  transport  is 
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(8) 


where  N  is  the  superficial  molar  flux  of  cations  through  the  porous 
region  and  is  calculated  by 


dc 


(9) 


where  is  the  cationic  transference  number.  The  appropriate  effec¬ 
tive  diffusivity,  Deff,  for  each  porous  layer  was  used.  The  boundary 
conditions  at  the  electrodes  (x  =  0,  L)  are  N  =  I/F.  Also,  electrolyte 
concentration  and  superficial  molar  flux  are  continuous  at  the 
porous-layer  interfaces. 


4.1.3.  The  intrinsic  transport  properties  of  the  electrolyte 

The  intrinsic  transport  properties  (D,  /c,  and  t°)  used  in  this 
model  are  obtained  directly  or  estimated  from  the  literature. 
Recently,  Stewart  and  Newman  [23]  measured  the  salt  diffu¬ 
sion  coefficient  for  LiPF6  in  1:1  (w:w)  EC/DEC  by  the  method 
of  restricted  diffusion  with  concentration  measured  by  ultravio¬ 
let/visible  absorption.  Their  fit  of  the  results  is 

D  =  2.582  x  10“5  exp  (-2.856C),  (10) 

where  D  is  in  units  cm2  s-1  and  c  is  in  units  mol  dm-3. 

For  conductivity  of  the  electrolyte,  we  used  a  concentration- 
dependent  shape  function  [18],  scaled  by  an  experimental  value 
of  the  conductivity  at  a  single  concentration,  namely  /cref  = 
7.8mScnrr1at  cref  =  1  moldm-3for  LiPF6  in  EC:DEC  [28].  The 
resulting  electrolyte  conductivity  used  in  this  model  is 

1.262(c/cref) 

1  +  0.2(c/cref)2  +  0.08(c/crej 

Note  that  the  small  non-zero  value  of  /c  at  zero  concentration  was 
used  in  the  model  to  ensure  numeric  stability  generally,  though 
in  this  particular  work  modeled  concentrations  never  approached 
zero. 

As  for  cationic  transference  number,  we  used  a  fixed  value,  = 
0.36,  in  accordance  with  the  earlier  work  employing  this  electrolyte 
[18]. 


+  0.014 


(11) 


4.2.  AC  impedance  experiment 


At  high  frequencies  the  concentration  in  an  AC  impedance 
experiment  is  effectively  constant,  which  effects  a  significant 
simplification  in  the  model.  In  addition,  one  can  neglect  the 
charge-transfer  impedance  at  the  electrode  surfaces  due  to  the  low 
impedance  of  double-layer  capacitance.  This  means  that  can  be 
directly  related  to  electrolyte  conductivity: 


R 


oo 


(12) 


where  Lj  and  Kj >eff  are  the  thickness  and  effective  conductivity  of 
layer  j.  The  summation  is  over  all  porous  layers  between  the  two 
electrodes.  A  is  the  area  of  the  electrodes  normal  to  the  current  flow. 


Fig.  5.  Illustration  of  possible  current  paths  in  a  cathode  film.  The  current  flows 
mostly  through  the  least-resistance  path  (arrows),  making  unreliable  the  measure¬ 
ment  of  Ke ff  without  detailed  knowledge  of  the  local  behavior  of  electronic  and  ionic 
paths. 


Eq.  (1)  is  used  to  calculate  the  tortuosity  of  the  separator  layers 
from  the  calculated/measured  /cj  eff  values  and  the  intrinsic  ionic 
conductivity. 

The  above  formula  does  not  work  well  for  electronically  con¬ 
ductive  cathode  films  as  used  in  our  experiments.  In  the  limit  of 
high  frequency  and  large  interfacial  area,  the  resistance  for  cur¬ 
rent  flow  between  ionic  and  electronic  pathways  in  the  cathode 
becomes  negligible.  Thus,  the  alternating  current  is  carried  by  both 
pathways,  ionic  and  electronic,  as  illustrated  in  Fig.  5.  The  figure 
shows  a  resistive  circuit  representing  the  cathode  film  between 
two  separator  layers.  The  current  will  follow  the  least-resistance 
path,  which  may  involve  multiple  hops  between  the  ionic  and  elec¬ 
tronic  modes  of  transport.  Therefore,  determination  of  /ceff  is  not 
straightforward. 

As  a  further  example,  consider  the  continuum  result  for  the  same 
one-dimensional  geometry  of  Fig.  5: 


Roo  —  Rsep  +  ^surf  + 
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(14) 

(15) 


where  o  is  the  electronic  conductivity  and  Rsurf  is  an  interfacial 
resistance,  which  will  be  explained  later  in  this  work  (Eq.  (17)). 
Functions  /c(x)  and  cr(x)  are  the  local  variations  in  ionic  and  elec¬ 
tronic  conductivities  respectively,  The  subscript  ‘avg’  is  used  to 
indicate  respective  average  quantities.  What  we  seek  is  /cav g,  which 
is  /ceff,  in  order  to  compute  the  tortuosity.  What  the  high-frequency 
AC  impedance  experiment  provides  by  Eq.  ( 13 )  is  (/c  +  cr)avg.  It  is  also 
possible  to  perform  dry  conductivity  measurements  on  the  cath¬ 
ode  film  to  obtain  <Tavg.  Elowever,  if  local  variations  in  ionic  and 
electronic  conductivities  are  negatively  correlated  (as  we  expect 
from  physical  arguments)  then  (/c  +  cr)avg  >  (/cavg  +  <ravg).  In  sum¬ 
mary,  because  the  local  correlations  of  conductivities  /c(x)  and  cr(x) 
are  unknown,  we  make  no  attempt  in  this  work  to  extract  /cavg  and 
hence  r  from  R ^  values.  We  nevertheless  believe  there  is  impor¬ 
tant  information  that  can  be  obtained  from  such  experiments  when 
combined  with  appropriate  modeling,  but  such  is  beyond  the  scope 
of  the  present  work. 

The  problem  just  outlined  is  not  present  for  separator-only 
experiments  where  a  ^  0.  Therefore,  we  used  AC-impedance  to  find 
the  tortuosity  of  the  separator  layers  and  to  serve  as  a  validation 
for  the  polarization  interrupt  experiment  on  the  same  system.  The 
polarization  interrupt  experiment  is  then  used  solely  to  determine 
the  tortuosity  in  the  active-material  films. 
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5.  Results  and  discussion 

5.1  Tortuosity  in  the  separator 

5.1.1.  AC  measurements 

Because  the  conductivity  measurements  were  done  at  high  AC- 
signal  frequency,  there  is  negligible  concentration  variation  in  the 
cell:  c  =  cre f.  Therefore,  we  take  k  =  /cref.  From  this  and  the  effective 
conductivity  (/ceff),  the  apparent  tortuosity  (rapp)  for  the  separator 
was  calculated  using  Eq.  (1).  This  experiment  was  repeated  using 
different  numbers  of  separators  and  the  apparent  tortuosity  of  the 
separators  was  calculated  in  each  instance.  Then,  the  results  for  the 
tortuosity  for  different  numbers  of  separators  were  compared  with 
those  obtained  from  the  polarization-interrupt  (Fig.  7). 


5.1.2.  Restricted- diffusion  measurements 

Fig.  6  shows  results  on  a  semi-log  plot  from  a  representative 
experiment  containing  three  separator  layers,  compared  to  the  cor¬ 
responding  model  results.  The  cell  was  polarized  for  2  min  at  a 
constant  current  density  (1.25  mAcnrr2)  and  then  the  current  was 
interrupted  for  10  min,  during  which  the  cell  potential  approached 
zero.  As  discussed  previously,  the  procedure  was  repeated  for  dif¬ 
ferent  current  densities  and  polarities.  The  tortuosity  in  the  model 
was  varied  to  best  reproduce  the  experimental  relaxation  time  con¬ 
stant  (slope  on  the  semi-log  plot).  The  best  tortuosity  was  taken  to 
be  the  one  that  minimized  the  sum  of  squared  error  during  the 
first  2  min  of  the  relaxation  period.  The  same  fitting  procedure 
was  repeated  for  different  current  densities  to  validate  and  aver¬ 
age  the  results.  The  model  curves  showed  good  agreement  with 
corresponding  experimental  curves,  as  illustrated  in  Fig.  6. 

Fig.  7  shows  the  apparent  tortuosity  (rapp)  as  a  function  of  num¬ 
ber  of  the  separators  (N),  as  obtained  from  AC-impedance  and 
polarization-interrupt  experiments.  One  can  think  of  the  verti¬ 
cal  axis  (rappN)  as  a  dimensionless  total  resistance  between  the 
electrodes.  The  first  thing  to  note  is  the  satisfactory  agreement 
between  the  AC-impedance  and  restricted-diffusion  experiments. 
Also  shown  in  the  figure  is  a  linear  fit  for  the  combined  experimental 
data  set.  The  slope  of  the  fit  represents  the  actual  tortuosity  of  the 
separator  layers,  calculated  here  for  Celgard  2400  to  be  r  =  3.15. 
However,  the  apparent  tortuosity  measured  for  one  separator  is 
tapp  =  6.12  which  is  in  good  agreement  with  rapp  =  5.98  obtained 
by  Patel  et  al.  [12].  It  is  perhaps  surprising  that  r  is  so  large  for 
this  material  given  that  the  pores  are  sometimes  represented  as 
being  straight  and  parallel.  Instead,  r  >  1  and  the  effective  transport 
properties  are  much  reduced  relative  to  that  ideal. 

5.1.3.  Interfacial  resistance 

Interestingly,  the  intercept  of  the  fit  in  Fig.  7  does  not  go  to  zero. 
This  suggests  that  there  was  an  additional  resistance  in  the  sys- 


Fig.6.  Model-experiment  comparison  for  galvanostatic  polarization  (t  <  120  s),  fol¬ 
lowed  by  interrupt  and  relaxation  (t  >  120  s)  in  a  cell  containing  three  separators. 


Fig.  7.  The  variation  of  the  apparent  tortuosity  rapp  with  the  number  of  separators. 
The  data  suggest  good  agreement  between  the  AC-impedance  and  polarization- 
interrupt  experiments.  The  slope  of  the  linear  fit  corresponds  to  the  true  tortuosity, 
while  the  non-zero  intercept  suggests  a  surface  resistance  at  the  lithium  electrodes. 


tern,  not  associated  with  the  separators  themselves,  but  instead 
associated  with  the  surfaces  of  the  lithium  electrodes.  This  inter¬ 
facial  resistance  could  be  due  to  the  solid  electrolyte  interphase 
(SEI)  or  some  kind  of  interfacial  spreading  resistance  or  both.  For 
instance,  active  redox  sites  on  the  electrode  surface  are  not  nec¬ 
essarily  aligned  with  the  separator  pores  and  spreading  resistance 
could  be  incurred.  Also  of  note  is  that  this  interfacial  resistance  has 
two  aspects:  ohmic  resistance  and  diffusional  resistance,  as  given  by 
the  fact  that  the  AC-impedance  and  polarization-interrupt  results 
have  essentially  the  same  intercept  in  Fig.  7.  Additional  work  is 
needed  to  determine  any  electrolyte-composition  dependence  of 
this  resistance. 

There  are  multiple  ways  one  could  take  into  account  this  addi¬ 
tional  interfacial  resistance.  It  is  not  accurate,  however,  to  model  it 
solely  as  an  electrical  resistance,  as  this  would  neglect  the  concen¬ 
tration  polarization  that  happens  as  well.  In  the  present  model  we 
included  the  interfacial  resistance  by  taking  it  to  be  a  thin  surface 
layer,  essentially  a  porous  separator,  placed  next  to  each  electrode. 
To  do  this  we  assumed  a  thickness  of  the  layer,  Lsurf  =  2.5  pan,  and  a 
porosity,  esurf  =  csep  =  0-37.  The  intercept  t*  of  Fig.  7  and  the  thick¬ 
ness  of  the  normal  separator  layer  Lsep  are  then  used  to  calculate  an 
effective  tortuosity  of  the  surface  layer: 


^  r*  ^sep^surf  /i  q 

rsurf  =  r  57 - ~ -  =  14.8. 

■^surf^sep 


(16) 


The  particular  values  used  for  Lsurf  and  esurf  are  not  important  as 
long  as  this  additional  surface  layer  has  the  appropriate  overall 
resistance  and  occupies  a  negligible  volume  fraction  of  the  cell, 
and  numerical  stability  can  be  maintained  in  the  model.  Because 
we  model  the  interfacial  resistance  as  a  distinct  layer  next  to 
each  electrode,  we  use  the  true  tortuosity  of  the  separators  (the 
slope  in  Fig.  7)  rather  than  the  apparent  tortuosity  to  model  the 
polarization-relaxation  experiments  of  the  cells  containing  cathode 
films.  Finally  we  note  that  in  the  case  of  the  AC  impedance  exper¬ 
iment  represented  in  Eq.  (13),  the  total  high-frequency  resistance 
associated  with  the  two  electrode  surface  layers  is 


Rsur[A  =  2LsurfTsurf  =  2.56  £2  cm2. 
KCsurf 


(17) 


5.2.  Tortuosity  in  the  active -material  films 


To  better  understand  the  effect  of  morphology  and  porosity  on 
the  tortuosity  of  cathode-material  films,  the  polarization-interrupt 
experiment  was  performed  on  cells  of  different  compositions  and 
porosities.  The  main  focus  in  this  paper  is  on  the  relationship 
between  porosity  and  tortuosity  for  LiFeP04  cathodes,  with  a  goal 
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Time  (s) 

Fig.  8.  Model-experiment  comparison  for  polarization  and  relaxation  in  the  LiFeP04 
film.  Two  model  curves  are  shown:  the  first  is  considered  the  best  fit  of  the  relaxation 
data.  The  second  is  shown  to  illustrate  the  sensitivity  of  the  model  result  to  tortuosity 
by  using  a  value  that  differs  from  the  best  fit  by  5%. 

to  develop  an  empirical  correlation  (like  the  Bruggeman  relation) 
between  these  two  variables.  However,  some  additional  cells  con¬ 
taining  LiCo02  active  material  were  tested  as  well. 

5.2.1.  Restricted-diffusion  measurements 

As  with  the  separator-only  cells,  the  detailed  model  is  used  to 
regress  the  tortuosity.  In  this  case  the  tortuosity  of  the  separator 
layers  (and  electrode  surface  layers)  are  fixed  and  the  tortuosity 
of  the  cathode  film  is  varied  to  best  fit  the  time  constant  of  the 
relaxation  curves. 

Fig.  8  shows  a  representative  semi-log  plot  of  model-experiment 
comparison  of  the  polarization-relaxation  experiment  for  a  LiFeP04 
film  (115  |jim  thick,  e  =  0.45).  The  cell  was  polarized  for  2  min  with 
a  constant  current  density  (1.25  mA  cm-2 )  and  then  the  curre  nt  was 
interrupted,  causing  the  cell  potential  to  approach  zero.  Shown  in 
the  figure  are  two  model  results.  The  first  curve  is  the  best  fit  of  r 
from  the  relaxation  time  constant  (slope  in  the  plot).  In  the  second 
model  curve  r  was  purposely  offset  from  the  optimum  by  5%  in 
order  to  illustrate  the  sensitivity  of  the  model  to  changes  in  r. 

A  statistical  analysis  was  performed  for  a  single  cell  to  deter¬ 
mine  a  characteristic  uncertainty  for  the  obtained  tortuosities.  The 
model  was  independently  fit  to  the  relaxation  results  for  five  dif¬ 
ferent  polarization  current  densities.  These  independent  samples 
yielded  a  standard  error  of  1.1%  in  tortuosity  at  the  90%  confidence 
level.  In  addition  to  this  measurement  uncertainty,  another  source 
of  error  is  the  cell-to-cell  variation  due  to  imperfect  fabrication  and 
natural  material  differences,  as  apparent  in  Fig.  9. 

5.2.2.  Tortuosity-porosity  correlation 

Accurate  cell  modeling  for  prediction  and  optimization  requires 
accurate  transport  properties.  To  aid  in  this  type  of  modeling  we 


Fig.  9.  The  effect  of  porosity  on  the  tortuosity  of  cathode  films.  The  tortuosity  is  well 
correlated  by  a  power-law  equation  (inset). 


present  Fig.  9,  which  summarizes  the  effect  of  porosity  on  the  tor¬ 
tuosity  of  the  tested  cathode  films,  with  the  focus  being  on  LiFeP04 
films.  A  power  function  is  fit  to  the  data  and  shown  in  the  figure. 
The  root-mean-square  fractional  deviation  of  experimental  points 
from  the  fit  is  2.7%.  The  maximum  absolute  fractional  deviation  for 
any  point  is  4.7%.  This  means  that  the  experimental  points  are  well 
described  by  the  correlation 

r  =  1.8c-0'53.  (18) 

This  means  a  =  1.53  in  Eq.  (4),  which  is  nearly  same  as  that  rec¬ 
ommended  by  Bruggeman  [14].  However,  the  tortuosity  obtained 
in  this  work  is  nearly  two  times  {y  =  1.8)  higher  than  that  pre¬ 
dicted  by  the  Bruggeman  relation.  This  suggests  that  the  effective 
transport  properties  can  be  well  predicted  by  a  scaled  Bruggeman 
relationship.  The  scaling  parameter  y  can  be  adjusted  for  differ¬ 
ent  constituents  and  particle  morphologies  of  the  active-material 
film.  Nevertheless,  in  Fig.  9  the  few  results  for  LiCo02 -material 
films  (indicated  by  triangles)  and  carbon-fiber-containing  LiFeP04- 
material  films  (indicated  by  squares)  show  good  agreement  with 
the  baseline  LiFeP04  results.  In  these  cases  either  the  active- 
material  differs  with  respect  to  particle  shapes  and  sizes,  or  the 
carbon  aspect  ratio  differs,  but  the  relative  mass  fractions  of  active 
material,  carbon  additive,  and  binder  are  the  same. 

Some  researchers  have  used  the  Bruggeman  exponent  as  an 
adjustable  parameter  in  full  battery  simulations  and  arrived  at  a 
relatively  large  value  [10,13].  However,  it  is  possible  that  the  large 
exponent  compensates  for  assuming  a  unity  scaling  factor  or  other 
possible  model  deficiencies.  Prior  work  by  our  group  [18]  used  the 
normal  Bruggeman  exponent  {a  =  1 .5);  however,  the  scaling  factor 
was  much  larger  {y  =  21)  than  obtained  here.  This  was  needed  to 
obtain  apparent  concentration-polarization  effects  observed  for  a 
thick  cell  under  high-rate  discharges.  However,  that  work  did  not 
use  a  concentration-dependent  liquid-phase  diffusivity.  Using  the 
diffusion  correlation  by  Stewart  and  Newman  [23],  in  which  D  can 
take  quite  small  values  at  large  concentrations,  better  enables  the 
model  to  reproduce  high-rate  discharges  for  thick  cathodes  with¬ 
out  resorting  to  such  a  large  scaling  factor.  In  addition,  the  earlier 
work  did  not  entail  a  direct  measure  of  liquid-phase  transport.  We 
therefore  believe  the  present  work  to  be  a  more  reliable  indication 
of  the  tortuosity  in  typical  porous  cathodes  for  lithium-ion  batter¬ 
ies.  The  Bruggeman  relationship  has  been  shown  to  be  accurate  for 
a  conductive  phase  mixed  with  monodisperse  insulated  spherical 
particles  [12,14,16,17].  TEM  images  of  LiFeP04  (Fig.  10)  suggest  a 
non-monodisperse  and  non-spherical  particle  morphology  for  the 
system.  The  addition  of  carbon  additives  (carbon  black,  graphite, 
and  carbon  fibers)  in  such  cells  further  distorts  the  pore  structure 
from  the  Bruggeman  ideal.  It  appears  that  the  Bruggeman  relation  is 
a  limiting  case  for  particulate  systems;  realistic  (non-monodisperse 
and  non-spherical)  systems  will  have  greater  tortuosity  than  given 
by  Bruggeman. 

5.2.3.  Model  validation  with  different  thicknesses 

To  optimize  batteries  for  HEV  and  PHEV,  it  is  also  desirable  to 
know  the  effect  of  different  resistances  with  the  thickness  of  the 
electrode  [29].  Our  present  work  helps  to  predict  the  liquid  phase 
resistance  with  different  thicknesses  of  the  electrode.  By  incor¬ 
porating  the  tortuosity  in  the  separator,  interfacial  resistance  on 
lithium  electrodes,  and  tortuosity  of  the  active-material  film,  the 
model  can  be  used  to  predict  the  diffusive  time  constant  for  dif¬ 
ferent  active-material  film  thickness.  The  results  shown  in  Fig.  9 
at  e  =  0.45  are  for  cathode  films  around  100  p,m  thick.  Fig.  11 
shows  results  for  a  thicker  (140  p,m,e  =  0.45)LiFePO4  film,  with  the 
model  operated  in  predictive  mode.  The  agreement  between  model 
prediction  and  experiment  is  good  in  terms  of  the  correct  relax¬ 
ation  slope  on  the  semi-log  plot.  This  supports  the  validity  of  our 
tortuosity-porosity  correlation  (Eq.  (18))  for  different  thicknesses. 
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Fig.  10.  TEM  image  of  LiFeP04  particles  suggesting  a  non-monodisperse  and  non- 
spherical  particle  distribution. 
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Fig.  11.  Transport  model  validation  from  a  polarization-interrupt  experiment  for  a 
thicker  LiFeP04  cathode  film  ( 140  p,m)  at  45%  porosity,  where  the  effective  transport 
properties  in  the  model  are  parametrized  from  experiments  on  thinner  (s»100  p,m) 
cells.  The  good  agreement  between  the  model  and  experiment  in  terms  of  relaxation 
decay  rate  suggests  the  tortuosity  correlation  given  in  the  text  is  accurate  for  different 
thicknesses. 

6.  Conclusion 

In  order  to  have  a  model  that  can  accurately  predict  cell  per¬ 
formance  under  varying  fabrication  parameters  and  operating 
conditions,  it  is  imperative  to  have  accurate  transport  properties 
in  the  model.  In  this  study,  experimental  and  modeling  methods  to 
determine  the  tortuosity  in  the  separator  and  the  active-material 
film  were  developed.  The  tortuosity  is  a  geometric  parameter  inde¬ 
pendent  of  the  particular  electrolyte  used;  its  value  is  needed  to 
predict  effective  conductivity  and  diffusivity  in  a  porous  mate¬ 
rial. 

AC  impedance  and  a  polarization-interrupt  (or  restricted  diffu¬ 
sion)  experiment  with  a  special  cell  geometry  were  used  to  measure 
directly  the  tortuosity  in  the  separator  and  the  active-material  films. 
Nevertheless,  and  particularly  in  the  case  of  the  active-material 
films,  interpretation  of  the  experiments  requires  a  detailed  math¬ 
ematical  model.  The  tortuosities  for  the  separator  obtained  from 
the  two  types  of  experiment  were  consistent  with  each  other.  This 
suggests  that  a  restricted-diffusion-type  experiment  can  be  used 
reliably  to  determine  the  tortuosity  of  a  porous  structure.  There¬ 
fore,  the  tortuosity  of  the  active-material  films  was  extracted  from 


the  model-experiment  comparison  of  the  polarization-interrupt 
experiments. 

Our  experiments  for  the  separator  uncovered  an  interfacial  resis¬ 
tance  associated  with  the  surface  of  the  lithium  electrodes  that 
prior  researchers  may  have  inadvertently  included  in  the  resistance 
assigned  to  the  separator,  hence  our  use  of  the  term  “apparent  tortu¬ 
osity.”  The  apparent  tortuosity  we  obtained  for  one  separator  layer 
was  found  to  be  in  good  agreement  with  values  reported  by  others. 
A  separate  accounting  for  the  interfacial  layer  leads  to  the  true  tor¬ 
tuosity  of  the  separator;  we  expect  this  value  to  be  more  appropriate 
in  cell  modeling. 

Also  interesting  is  that  our  experiments  suggest  this  interfacial 
layer  restricts  salt  diffusion  with  the  same  tortuosity  as  it  does  ionic 
current,  at  least  for  the  particular  electrolyte  used  here.  The  surface 
resistance  can  be  modeled  as  an  additional  porous  layer  adjacent 
to  each  lithium  metal  electrode.  Our  prior  experience  suggests  con¬ 
centration  depletion  can  occur  in  lithium  batteries  cycled  at  high 
rates  («  5  C)  and  thus  have  a  significant  effect  on  cell  voltage.  There¬ 
fore  it  is  important  to  correctly  account  for  diffusive  resistances  in 
separators  and  surface  films  as  well  as  in  the  porous  electrodes.  In 
addition,  one  must  correctly  account  for  the  concentration  depen¬ 
dence  of  the  liquid  diffusivity. 

A  correlation  between  tortuosity  and  porosity  of  active-material 
films  was  obtained.  The  oft-used  Bruggeman  exponent  {a  =  1.5) 
seems  to  predict  the  porosity  dependence  of  tortuosity  quite  well. 
However,  the  tortuosities  of  the  porous  structures  studied  in  this 
work,  both  the  separator  and  active-material  films,  are  approxi¬ 
mately  two  times  larger  than  those  predicted  by  the  Bruggeman 
relation.  In  particular  for  the  separator  this  confirms  that  the  trans¬ 
port  paths  are  not  straight  parallel  channels. 

The  experimental  techniques  and  models  developed  in  this  work 
can  be  used  to  determine  the  transport  resistances  for  cathode  films 
with  different  porosities  and  thicknesses.  This  can  help  in  optimiz¬ 
ing  batteries  for  different  applications.  Once  the  factors  that  affect 
the  different  transport  resistances  are  quantified,  fabrication  strate¬ 
gies  can  be  developed  to  mitigate  the  dominant  resistances.  For 
instance,  lower  porosity  is  desired  for  lower  electronic  resistance 
[19].  On  the  other  hand,  higher  porosity  or  liquid-volume  frac¬ 
tion  is  desired  for  lower  liquid-phase  resistance.  Therefore,  there 
is  a  trade-off  between  electronic  and  liquid-phase  resistance  with 
porosity  [29].  Once  the  electronic  resistance  is  known,  our  model 
and  experimental  methods  can  be  used  to  predict  such  an  optimum 
point  for  cell  fabrication  for  a  particular  application. 
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